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ABSTRACT 

Seasonal variations of firn density in ice-sheet firn layers have been attributed to 
variations in deposition processes or other processes within the upper firn. A recent high- 
resolution (mm scale) density profile, measured along a 1 8 1 m core from Antarctica, 
showed small-scale density variations with a clear seasonal cycle that apparently was not- 
related to seasonal variations in deposition or known near-surface processes (Gerland and 
others 1999). A recent model of surface elevation changes (Zwally and Li, submitted) 
produced a seasonal variation in firn densification, and explained the seasonal surface 
elevation changes observed by satellite radar altimeters. In this study, we apply our 1-D 
time-dependent numerical model of firn densification that includes a temperature- 
dependent formulation of firn densification based on laboratory measurements of grain 
growth. The model is driven by a steady-state seasonal surface temperature and a 
constant accumulation rate appropriate for the measured Antarctic ice core. The modeled 
seasonal variations in fim density show that the layers of snow deposited during spring to 
mid-summer with the highest temperature history compress to the highest density, and 
the layers deposited during later summer to autumn with the lowest temperature history 
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compress to the lowest density. The initial amplitude of the seasonal difference of about 
0.13 reduces to about 0.09 in five years and asymptotically to 0.92 at depth, which is 
consistent with the core measurements. 


INTRODUCTION 


Seasonal variations of fim density have been observed in numerous field snow-pit and 
fim-core studies. These variations are more pronounced in the top several meters fim, 
and generally becoming less significant with increasing depth. Interpretation of the 
causes of these density variations during fim diagenesis is usually complicated by the 
variety of changing surface weather conditions. At Plateau station in Antarctica where 
there is no summer melting, the summer section of each annual fim layer generally 
consisted of a thin, hard, fine-grained layer, while the winter section was softer, coarse- 
grained, and more homogeneous than the summer section (Koemer, 1971). However, 
snow-stratigraphy studies at other polar locations generally showed that summer layers 
are normally coarser-grained and have lower density and hardness values (e.g. Gow 
1961; Benson 1962). 

Using improved technology, a very detailed density profile was continuously measured in 
3 mm depth intervals along a 181 m core (Gerland and others, 1999). This fim/ice core 
was from the summit of Berkner Island (79° 3651 S, 45° 4328 ), Antarctica, which has 
an annual mean temperature of -27 °C and insignificant summer melting. Such high- 
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resolution measurement provides an excellent and important example for investigating 
the small-scale characteristics of density evolution during dry-fim densification, which 
previously has not been satisfactorily explained. In particular, Gerland and others (1999) 
noted that the measured seasonal cycle, which shows higher densities of summer layers 
found in the depth below 25 m, was not consistent with observations of higher winter 
densities in some locations. 

Meanwhile, recent analysis of satellite radar altimetry measurements of ice sheet surface 
elevations revealed seasonal variations in dry firn zones. The observed winter surface is 
generally higher than the summer surface, with the minimum elevation typically 
occurring before or near the minimum in summer temperature (Zwally and Li 
Submitted). This cycle has not been satisfactorily explained either on the basis of 
seasonal variations in the effective depth of the radar-measurement redundant or by 
seasonal variations in snow accumulation at the surface. Also, the temperature sensitivity 
of previous densification models was insufficient to produce a seasonal elevation cycle of 
sufficient amplitude (e.g. Yi and others 1997). 

Our model of surface elevation changes (Zwally and Li, submitted) explained the 
seasonal surface elevation cycle observed by satellite radar altimeters and produced a 
significant seasonal variation in firn density profiles. In this study, we apply our 1-D 
time-dependent numerical model to simulate detailed variations in firn density measured 
in the Antarctic core by Gerland and others (1999). A principal feature of the model is 
the temperature-dependent formulation of firn densification that is based on laboratory 
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measurements of grain growth. We use a steady-state surface temperature with seasonal 
amplitude estimated from a nearby AWS and a constant accumulation rate to drive the 
model. 

THE MODEL 

Essential to modeling of the density-depth profile is the rate equation for fim 
densification. The empirical relation (equation (1)) given by Herron and Langway (1980) 
is used in the model for calculating the densification rate. 

dp/dt =kA a (pi-p) (1) 

where k is a “rate constant” that is solely dependent on fim temperature (T), A is the 
mean accumulation rate and exponent a is a constant depending on the stage of 
densification. Equation (1) has been based on the suggestion that in the densification of 
fim, the proportional change in air space is linearly related to the change in stress due to 
the weight (p) of overlying snow (Robin 1958). In equation (1) this weight change is 
represented by accumulation rate (A) which is normally a function of time (t) i.e. dp/dt = 
A(t). 

The temperature dependence of k follows the Arrhenius-type relation shown by 

k = k 0 exp(- E/RT) (2) 
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where R is the gas constant. Both ko and the activation energy E have often been used as 
constants independent of temperature. However, field and laboratory experiments for 
grain growth and ice creep indicate that the activation energy E for each process is 
actually a function of temperature (Jacka and Li 1994). Zwally and Li (submitted) 
incorporated this feature into the model by using the values of E (T) obtained from a best 
fit curve through the data of activation energy for grain growth. Corresponding ko(T) 
values are used to obtain k(T) consistent with the measurements. 

The link between the changes in depth and density is given by equation (3). A fim layer 
is considered with initial thickness determined by the precipitated snow with surface 
density, p 0 in a time span, dt. Equation (1) is applied to each fim layer to calculate the 
densification rate during dt and thus the density. Assuming the density within each layer 
is reasonably uniform, for the layer at a depth i and time j the corresponding change (Vjj ) 
of the layer thickness (hjj) caused by the density change, is given by 

vq = hjj/pij (dp/dt) ; j (3) 

The burial depth (z) at any given layer can then be estimated by the sum of the remaining 
thicknesses of all fim layers at the depth above. 

Variations of temperature with time in each fim layer is computed based on the standard 
heat transfer equation (Paterson 1994) 
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dT 

dt 



( 4 ) 


where f is the thermal diffusivity depending on fim density. We have also used the 
analytical solution of equation (4) assuming the surface air temperature varies in regular 
sine function of time with the amplitude T s (Paterson 1994). The temperature and the 
density for a fim layer were calculated at the same time. 


Equation (1) only takes into account the effects of stress and temperature. The usage of 
the equation (1) as the constitutive relation to describe the densification process for near 
surface fim (cm) is rather simplified. Previous studies show that besides other surface 
processes such wind, solar radiation, temperature gradient (dT/dz) may also accelerate 
the densification and grain growth in the near surface fim (Alley 1988; Colbeck 1989; 
Gray and Morland 1995). However these physical processes are more significant in low 
accumulation area because the intensity of these effects also depends on the duration 
(determined by accumulation rate) for which the surface fim is subjected to the highly 
variable surface weather conditions (Nishimura and Maeno 1984). In high accumulation 
area such as east side of Law Dome Antarctica, these enhanced rates are greatly reduced 
and the depth hoar layers are hardly developed (Xie and others 1989; Li and others 1990). 
Therefore our model may be more feasible in the dry snow area with higher accumulation 
rates where the snow stratigraphy profile is relatively uniform. 

RESULTS AND DISCUSSION 
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The model is driven by the surface air temperature with a regular seasonal cycle and by a 
constant accumulation rate. Time steps are 30-days. AWS records from two nearby 
stations AGO-A84 (84.36°S, 23.86°W) and Limbert (75.42°S, 59.95°W), respectively 
located at higher and lower latitude than the drilling site, show the maximum variation of 
monthly mean surface temperature is approximately 3 1 0 (Linda and others 1996, 1997). 
Therefore, we use 15.5° as the approximate amplitude (T s ) of the surface temperature 
with annual mean -27 °C to derive the fim temperature for the drilling site. We also use 
annual mean accumulation rate 174 kg m' 2 a' 1 from the drilling site with surface snow 
density 300 kg m' 3 (Gerland and others 1999). Herron and Langway (1980) empirically 
derived a value of a in equation (1) of approximately 1 for the density below 550 kg m' 3 , 
and we use this value as an approximation for the entire density range. 

The modeled depth-density profile is compared to the field observation in Figure 1 . The 
shaded area shows the envelope of modeled density variations with depth. Two dashed 
lines with markers in Figure 1 indicate the range of observed seasonal density variations 
derived from the mean density profile given by Gerland and others (1999), and a 
modified curve depth-density variation range for the site (see below). In their core 
analysis, Gerland and others (1999) found a reverse change of the density variation 
between 25-50 m along the general decreasing trend of the density variation with 
increasing depth (density), forming a second peak of the density variation around the 
depth of 50 m. In this analysis, we consider that this second peak of the density variation 
is abnormal, which is possibly due to the interannual changes of the surface weather 
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conditions. Changes in physical properties of ideal dry' firn during densification more 
likely follow a continuous trend. By following the general trend of Gerland and others 
(1999, Fig. 4), a modified curve showing the range of seasonal variation of measured 
density with depth for the site is given in Figure 2 as indicated. 

As shown by Figure 1 , our modeled density depth-profile follows the observation with 
similar range of density variation that declines with increasing density or depth. 

However, it is evident from Figure 2 that the modeled density variation does not match 
the observed curve perfectly. In particular, the modeled densification rate exceeds the 
observed rate before the density reaches approximately 500 kg m' 3 , although the modeled 
density appears to be within the range of observed density variations. Besides any 
imperfections in the model, it should be noted that the observed density variations may 
bias to its steady state profile due to the possible interannual changes of surface weather 
conditions. 

After the detailed examination of the density data, Gerland and others (1999) indicate that 
the variations in firn density shown by Figure 1 are actually seasonal. The magnitude of 
these seasonal variations decreases continuously with depth. This feature of seasonal 
variation in firn density has been well modeled in our results. Figure 3 presents in larger 
scale the modeled density variations for the first 10 m shown in Figure 1. The magnitude 
of the variation decreases with depth as indicated by Figures 1, 2, and 3. This 
characteristic of the seasonal variation in fim density also agrees well with the field 
findings at Plateau station in Antarctica (Koemer, 1971). 
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To examine the specific seasonality in the density variation, we follow the density 
evolution for each fim layer deposited at 30-day time intervals within a one-year surface 
temperature cycle, as shown by Figure 4. The densification history of three selected fim 
layers together with its corresponding temperature history following deposition is shown 
in Figure 5. As shown by Figure 4, fim layer zero (CO) is deposited in early spring 
approximately three months prior to the maximum summer temperature. This snow 
begins to densify at a higher rate as the temperature increases following deposition. Fim 
layer 5 (C5) is deposited in autumn about two months after the maximum summer 
temperature. Layer C5 densifies a little while the temperature is in mid-cycle after the 
deposition, as indicated by the small step of the density curve. However, with the fall in 
temperature toward winter, the densification is essentially quenched until the next warm 
season starts in spring. The second step in density of C5 is smaller than the first step for 
layer CO deposited in the warm season, because C5 is insulated some by the winter 
deposited snow. This autumn snow obtained the least density of any of the 30-day layers. 
Snow deposited in winter obtains a density between these two extremes as shown by the 
third diagram for layer Cl 1. 

The longer period histories of temperature for fim layers CO and C5 in Figure 6 show that 
after subsequent seasons the temperature histories of the maximum and minimum density 
layers are not much different. Clearly, the temperature history of a layer during the first 
seasonal cycle is the primary determinant of the initial densification, which then remains 
correspondingly high or low as densification continues. The densification history for the 
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12 seasonal layers in Figure 7 shows that layers CO to C3 and layers Cl 1 to C12 that are 
deposited when temperature is increasing (i.e. dTs/dt >0) obtain higher densities than 
layers C4 to C9 that are deposited when the temperature is decreasing (dTs/dt <0). 
Therefore, it is the spring through early summer layers that obtain the highest densities 
and the autumn through early winter layers that obtain the lowest densities. 

CONCLUSION 

The high resolution density profile of 181m core retrieved from the summit region of 
Berkner island Antarctica provides an excellent example of fim densification, especially 
the seasonal scale density variations detailed to seasonal scale. We believe that the 
model presented here provides a good description of the densification process that leads 
to the observed seasonal cycle in fim. The surface temperature cycle is sufficient to 
generate the seasonal variations in density profile that is observed in dry polar fim. The 
maximum densities occur in the firn layers deposited during late spring to early summer 
and minimum densities in those deposited in autumn to early winter. 
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FIGURE CAPTIONS 

Figure 1 . Comparison of modeled and observed depth-density profiles for Summit of 
Berkner island, Antarctica. Dark area shows the envelope of modeled density 
variation. Dashed lines with markers indicate the range of observed density 
variation modified from Gerland and others (1999). 

Figure 2. Modeled (solid line) and observed (broken line) range of density variation 

versus depth for Summit of Berkner island, Antarctica. Data shown by the broken 
line are modified from Gerland and others (1999). 

Figure 3. Seasonal variations in fim density shown by modeled depth-density profile for 
Summit of Berkner island, Antarctica. 

Figure 4. One year surface air temperature (30 day mean) cycle in each 30 day time 

span. The figure beside each marker indicates the number of the 30 day from the 
beginning of the year corresponding to the time of deposition of the modeled 
layers. 

Figure 5. Densification history of three selected fim layers together with the temperature 
history experienced by the corresponding fim layer (cf. numbers in Figure 4). 
Late-winter layer Cl 1 and early-spring layer CO obtain higher densities than late 
autumn layer C5 (lowest density). 

Figure 6. Modeled temperature profiles showing the decrease of amplitude of fim 

temperature experienced by layers (CO dashed line and C5 solid line) while the 


13 


Submitted to Annals of Glaciology. 34A1 19 


fim layers travel downwards as implied by the changes with time. After the first 
year the temperature histories are similar. 

Figure 7. Sequence of the density history for each snow layer indicated in Figure 4 
showing the seasonal variations of the density in each seasonal fim layer. 
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